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ABSTRACT 

We test whether there is a relation between the observed tidal interactions and Seyfert activity 
by imaging in HI twenty inactive galaxies at the same spatial resolution and detection threshold as 
the Seyfert sample. This control sample of inactive galaxies were closely matched in Hubble type, 
range in size and inclination, and have roughly comparable galaxy optical luminosity to the Seyfert 
galaxies. We find that only ~15% of the galaxies in our control sample are disturbed in HI, whereas 
the remaining ~85% show no disturbances whatsoever in HI. Even at a spatial resolution of ~10 kpc, 
none of the latter galaxie s show appreciab le HI disturbances reminiscent of tidal features. 



In a companion paper (Kuo et al. , 2008), we report results from the first systematic imaging survey 
of Seyfert galaxies in atomic hydrogen (HI) gas. We find that only ~28% of the eighteen Seyfert 
galaxies in that sample are visibly disturbed in optical starlight. By contrast, ~94% of the same 
Seyfert galaxies are disturbed spatially and usually also kinematically in HI gas on galactic scales of 
> 20 kpc. In at least ~67% and up to perhaps ~94% of cases, the observed disturbances can be traced 
to tidal interactions with neighboring galaxies detected also in HI. The dramatic contrast between the 
observed prevalence of HI disturbances in the Seyfert and control samples implicates tidal interactions 
in initiating events that lead to luminous Seyfert activity in a large fraction of local disk galaxies. 
Subject headings: galaxies: active — galaxies: Seyfert — galaxies: interactions — galaxies: structure 
— galaxies: ISM 



1. INTRODUCTION 

Active Galactic Nuclei (AGNs) are believed to be man- 
ifestations of the vigorous accretion of gas onto a central 
SuperMassive Black Hole (SMBH). The origin or source 
of this gas, and how it is brought into the sphere of infiu- 
ence of the SMBH, is one of the longstanding unresolved 
pro blems in AGN r esearch (for a recent brief review, see, 
e.g., [Martini , 2004). Only a very small fraction (roughly 
1%) of galaxies m the local universe exhibit luminous 
AGNs, although an appreciable fraction may exhibit non- 
stellar nuclear activity at some level. 

One of the most popular hypothesis for triggering lumi- 
nous AGNs is gravitational interactions between galax- 
ies. Observational studies of local AGNs have mostly fo- 
cussed on optically-selected samples of Seyfert galaxies. 
Seyfert galaxies are preferentially found among early-tpe 
spiral galaxies, and exhibit very broad emission lines 
(from their nuclei) compared with their inactive coun- 
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terparts. Although most Seyfert galaxies do not ap- 
pear to be visibly disturbed let alone interacting in op- 
tical starlight, some studies report that Seyfert galaxies 
more frequently possess projected or genuine neighboring 
galaxies (within a given angular separation) co mpared 
with matched samples of inactive galax i es (e. g., Stauf- 
ferl [19821 IDaha ri, 1984; Rafane lh etHj [19951 [Pultzm- 



Hacyan et al.[ 1999; Ko ulouridis et al.[ 20061). This may 
constitute indirect evidence for more nrequent interac- 
tions between Seyfert and their neighboring galaxies than 
their visual appearences would suggest. On the other 
hand, there also are studies that refute reports of any 
differences between the frequency of Seyfert and inac- 
tive galaxi e s with projected neighboring galaxies (e . g . 



MacKentyl [19891 IFuentes- W illiams Stocke 



mm de 



R obertis, Yee, Hayhoer [T998 ; Sch mitt[ [20QT] ) 

Studies of nearby galaxies have demonstrated that 
atomic hydrogen (HI) gas can r eveal tidal features no t 
visible in optical starlight (e.g., jYun, Ho, fc Lo 1994). 
Because the HI disk of normal spiral galaxies usually ex- 
tends nearly twice as far out as the stellar optical disk, 
the outer regions of the HI disk is more loosely grav- 
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itationally bound. The HI disk is therefore more sus- 
ceptible to external gravitational perturbations than the 
stellar optical disk, and when perturbed the outskirts 
also take longer to dynamically relax. This makes HI 
the most sensitive and enduring tracer known of gravita- 
tional interactions between galaxies. To directly address 
whether Seyfert galaxies are involved in galaxy-galaxy in- 
teractions, we have therefore, for the first time, imaged 
uniformly a relatively large sample of Seyfert galaxies in 
HI gas. The full resu lts of that s tudy are presented in a 
companion paper by Kuo et al.l (12008). Here, we com- 
pare the results of Kuo et al. (2008| with those obtained 
for a matched sample of inactive galaxies; i.e., a control 
sample. 

To select a control sample, we apply the lessons learnt 
from optical studies that address whether galaxy-galaxy 
interactions are responsible for trigger ing Seyfert activ- 



ity (see, e.g., discussions in gFuente s- Willia ms fc Stocke 
T988llde Robertis, Yee, Hayhoel 1998; D ultzin-Hacyai^ 



et al 



1999). A number of the especially earlier studies 



select as their control sample the closest (in projected 
separation) inactive gal axies with sizes comparable to 



Rafanelli et al. 1995). As pointed by pTultzin-Hacyan et 



the Seyfe rt galaxy (e .g., |Dahari[ |1984| lMacKenty [ 1 1989 

IDul 



al.|(|1999), if Seyfert galaxies lie at or close to the center 
of a L 



ot aTocal galaxy density enhancement, a control sample 
selected in this manner may then systematically lie in a 
region o f lower galaxy den sity. Consistent with this pos- 
sibility, MacKenty (1989) find that a larger fraction of 
Seyferts have projected neighboring galaxies than com- 
parably sized inactive galaxies lying in the relatively close 
vicinity, but not with those selected at random across the 
sky. 

The more recent optical studies select at random in- 
active galaxies with the same morphological (i.e., Hub- 
ble) type and redshift, and sometimes also size and ab- 
solute magnitude, as the Seyfert galaxies. This match- 



ing is made either o n a on e-to-one basis (e.g., de Rober- 

or over the range of prop 



1998^ 



tis, Hayhoe fc Yee 

erties observed (e.g., Dultzin-Hacyan et al.| 1999). As 



poi nted out by both Fuentes- Willia ms fc Stockej ( 1988 ) 
and Dultzin-Hacyan et al.| ([1999| , matching in abi 
magnitudes may systematically bias the control sample 
to higher luminosities because of the non-negligible con- 
tribution from the AGN. An exact match to the absolute 
magnitude of the Seyfert galaxy minus its AGN contri- 
bution is usually not possible as such measurements are 
not generally avail able, but is desirable to the degree 
possible (e.g., as in de Robertis, Yee, fc HayhQe{ p998). 
By minimizing the differences in galaxy properties be- 
tween the Seyfert and control samples, it is hoped that a 
meaningful comparison of their environments will be pos- 
sible, alt hough this is not obviously the case . Using this 
method, de Robertis, Yee, fc Hayhoe (1998) find no sig- 
nificant excess of projected neighboring galaxies between 
their Seyfert and control samples, and at best only a 
marginal difference between their light asymmetries. On 
the other hand, Dultzin-Hacyan et al. (e.g., 1999) find 
an excess of Seyfert 2, but not Seyfert 1, with projected 
neighboring galaxies compared with their co ntrol sample, 
as wa s later confirmed spectroscopically by |Koulouridis| 
et al.l ([2006). 



In our HI imaging study, we select a control sample 
that is matched on a one-to-one basis in Hubble type. 



and which range from lower to approximately the same 
absolute magnitude as, the sample of active galaxies se- 
lected by ~" 



Kuo et al. 



(2008). We also matched over the 
observed range of optical sizes and inclinations. No vi- 
sual inspection of the control sample was made to see 
whether they are visibly disturbed, but we are aware that 
the selection criteria used may select against strongly 
disturbed galaxies not easily classified in Hubble type. 
Indeed, using much the same criteria to select their con- 
trol sample, de Robertis, Hayhoe fc Yee ('1 998) find that 
a higher proportion of Seyfert galaxies are involved in 
late-stage mergers than their control sample, although 
a similar fraction of the control sample displays signif- 
icant light asymmetries that could be evidence for re- 
cent interactions. Note tha t only a sm all fraction of the 
active galaxies studied by Kuo et al.l (12008) are visibly 
disturbed in the optical, but not so strongly that they 
cannot be classified in Hubble type. Nevertheless, as we 
shall see, the fractio n of optically-di sturbed galaxies in 
the Seyfert sample of Kuo et al. (2008*) is higher than that 



in our control sample. Because our control sample is ran- 
domly selected from the field, their environments should 
be representative (or at least not biased in any peculiar 
way against those) of inactive galaxies with similar Hub- 
ble types, absolute magnitudes, and sizes as the Seyfert 
sample of Kuo et a l. (2008). 

In §2 we describe how we selected our control sam- 
ple of inactive galaxies, and in §3 our HI observations of 
these galaxies and the data reduction. In §4 we present 
the results, and in §5 compare thes e results with thos e 
for the Seyfert sample studied by Kuo et al. (2008). 
In §6, we provide a concise summary of our princi- 
ple findings and their implications. Unless otherwise 
specified, we assume throughout a Hubble c onstant of 
Hp = 67 km s~^ Mpc~^ and ^7 = 1 as used by Kuo et al. 



(2008). 



2. SAMPLE SELECTION 



As described more fully in Kuo et al. (2008), the par- 
ent sample of active galaxies comprise all twenty- seven 
disk galaxies listed in the Veron-Ce tty fc Veron (1998, ), 
plus another from the Veron-Cetty fc Veron ( 200^, cat- 
alog at redshift s of 0.015 < z < 0.017 and with ab- 
solute 5-band magnitudes of —19 > Mb > —23 (for 
Ho = 50 km s~^ Mpc~^, as used in that catalog) in the 
northern hemisphere {S > 0°). Those galaxies were ob- 
served at an angular resolution of ~60'^ which at their 
distances of 66-75 Mpc corresponds to a spatial resolu- 
tion of ~20 kpc. The latter is roughly comparable with 
the optical disk sizes of the active galaxies, and hence 
the HI observations were tailored towards the detection 
of structures extended on galactic scales characteristic of 
tidal features. 

Our objective here is to conduct an equivalent HI imag- 
ing study of a comparable number of inactive galax- 
ies matched as closely as is reasonably possible to the 
active galaxy sample. We selected our control sample 



from th e CfA Redshift Survey (|Huchra et al.] [ 19831 p|995. 
see also [http: //cfa- www. harvard. edu^huchra/zcat/ ) ac- 



cording to the following criteria. 

1. Redshift: to reduce the required observing time 
(which amounted to an integration time of 2 hrs 
and an observing time of ~2.5 hrs per object in our 
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active galaxy sample), we selected inactive galax- 
ies at half the reshift range of our active galaxy 
sample; i.e., at 0.0075 < z < 0.0085 (distance of 
33.4-37.8 Mpc). A total of 166 galaxies in the CfA 
redshift catalog met this criterion. 

Individual samples selected in optical studies are 
located over a much broader range of redshifts, 
making one-to-one (or overall) matching in red- 
shift (range) important for avoiding difficulties 
when looking for an excess of projected neighboring 
galaxies between two samples. Here, our objective 
is to search for large-scale HI disturbances in our 
target galaxies, and determine whether these dis- 
turbances are produced by tidal interactions with 
neighboring galaxies. In such a case, matching in 
redshift is not necessary so long as we observe both 
samples at the same spatial resolution and detec- 
tion threshold (see §3), and provided there is no 
significant Cosmological evolution between the two 
redshifts. Of course, to look for interacting neigh- 
boring galaxies (but not necessarily to detect tidal 
features in the target galaxies), in both samples the 
field of view has to be sufficiently large to encom- 
pass any such galaxies. 

2. Optical morphology and luminosity: from the par- 
ent sample of 166 inactive galaxies, we matched 
on a one-to-one basis the Hubble type and optical 
5-band luminosity of each active galaxy to one or 
more inactive galaxies. Both the optical morphol- 
ogy and luminosity were taken fr om Leda, now re- 
named t h e Hyperleda database (|http:/ /leda.univ-| 



lyonl.fr; Paturel et al. , 2003). The morphology 



is matched to an accuracy of one Hubble subtype 
(e.g., for an Sa active galaxy, we matched to SOa- 
Sab inactive galaxies) given the inherent fuzziness 
in morphological classifications. Because the ab- 
solute magnitudes of the active galaxies include a 
poorly determined contribution from their bright 
nuclei, we allowed the match in absolute magni- 
tude to be up to one magnitude larger (i.e., optical 
5-band luminosity up to a factor of 2.5 lower) than 
that of the corresponding active galaxy. 

3. Optical size and inclination: based once again on 
the Hyperleda database, we also loosely matched in 
optical size by restricting the galaxy diameters to 
the range ~0^8-4^3 (with the majority in the range 

as measured at the 25 mag arcsec"^ 
isophotal level in the 5-band. This range in angu- 
lar diameters is roughly comparable with the corre- 
sponding range in physical diameters of the active 
galaxy sample. Because none of the active galaxies 
are close to face on (i.e., inclinations close to 0°), 
we also restricted the inclinations of the inactive 
galaxies to angles > 26° with respect to the plane 
of the sky. 

4. Radio continuum: to avoid the effect in our dy- 
namic range of our maps, we discarded those galax- 
ies that lie near bright radio continuum sources de- 
tected in the NRAS VLA Sky Survey (NVSS; Con- 
don et al. 1998) which would appear in the field 
of view of our HI observations. Here, we set a flux 



density cut at 21 cm of 1000 mJy. The detected 
flux densities of both sample and its companions 
are in the order of 100 mJy. This cut will not cause 
any bias when comparing results of the seyfert and 
control sample. 

In this way, we selected twenty-seven inactive galaxies 
that are closely matched in Hubble type and to the degree 
possible absolute magnitude, as well as range in size and 
in clination, t o the sample of twenty-eight active galaxies 
in|Kuo et al. (12008). The basic properties of this matched 
)le^( 



control sample^are listed in Table 1. Note that some of 
the values originally obtained from the Leda database 
have since been updated in the Hyperleda database; we 
list the latest, and presumably most reliable, values from 
the Hyperleda database. Our selection is entirely blind 
to the HI content of these galaxies. 

In Figure 1, we plot the number distribution in mor- 
phological types, absolute magnitudes, optical si zes, and 
inclinations of both the active galaxy sample in Kuo et 
al. (2008) and our control sample. The number distribu- 
tion in morphological types (top panel) shows no marked 
differences between the two samples within a precision of 
a Hubble subtype; a Kolmogorov-Smirnov (K-S) test ap- 
plied to the two samples gives a K-S statistic D = 0.18 
(converging to for identical distributions) and prob- 
ability Pr = 0.74 (1 for identical distributions). As in- 
tended, the absolute magnitudes (second panel from top) 
of the control sample are skewed towards larger values 
(i.e., lower luminosities) than the active sample. The 
distribution in optical sizes also is skewed towards some- 
what smaller dimensions for the control sample, but the 
majority in both samples span the same range. The con- 
trol sample has more moderately inclined (30°-45°) than 
high inclined (75°-90°) galaxies compared to the active 
sample, but once again the majority in both samples span 
the same range. The relatively small differences seen be- 
tween the distribution of physical properties between the 
active and control samples should not have a significant 
effect on the detecta bility of HI tidal features in the two 
samples; as shown by|Kuo et al.|(2008), the vast majority 
of the active galaxies, and nearly all those classified as 
Seyferts, exhibit HI disturbances usually in the form of 
tidal features irrespective of their Hubble type, absolute 
magnitude, optical size, or inclination within the selected 
range. 

3. OBSERVATIONS AND DATA REDUCTION 

We observed the control sample with the Very Large 
Array (VLA) of the National Radio Astronomical Obser- 
vatory (NRAO). Like the active galaxy sample studied by 



Kuo et al.| ( |2008|) , we used the most compact configura- 
tion of the VLA (the D-array). NGC 5375 had previously 
been adequately (for the purpose of this experiment) ob- 
served in this configuration on 1996 August 15, and so we 
simply retrieved the data for this galaxy from the VLA 
archive. The remaining galaxies were observed on 2003 
February 19, 21, and 24, apart from NGC 5355 which 
was observed on 2004 July 14. Other relevant details of 
the observations, including the flux, secondary (for am- 
plitude and phase calibration) , and bandpass calibrators 
used for each target object, are summarized in Table 2. 

The cor relator was confi gured in the same manner de- 
scribed in Kuo et al. (2008), recording signals in orthogo- 
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nal circular polarizations in sixty-four channels spanning 
a bandwidth of 6.25 MHz. The corresponding channel 
separation in velocity is ~21.2 km s~^ (the actual veloc- 
ity resolution is ~25.2 km s~^), and altogether span a ve- 
locity range of ~1350 km s~^. For the assumed Cosmol- 
ogy, this velocity range corresponds to a redshift interval 
of Az ~ 0.0045 (distance range of ~20 Mpc), compared 
with the redshift interval of Az ^ 0.001 for the control 
sample. The central channel is set to the systemic he- 
liocentric velocity of the HI line if previously measured 
(usually with a single-dish telescope), or if not then at 
its reported optical redshift. 

At 21 cm, the primary beam of the VLA is ~32' at 
full-width half-maximum, corresponding to a diameter 
of ~330 kpc at the distance of our control sample. Our 
field of view for the control sample is, of course, only 
half the linear diameter of our field of view for the active 
galaxy samp le, which is located at twice the distance. As 
described in Kuo et al. (2008), the projected separation 
between the Seyfert and their interacting neighboring 
galaxies is mostly (for ~85% of the sample) < 100 kpc. 
Our field of view is therefore sufficiently large to en- 
compass most of the interacting neighboring g alaxies de- 
tected around the Seyfert galaxies studied by |Kuo et al. 



(2008) even if the latter was placed at the distance of our 
control sample. 

We targeted the sa me detection thres hold at the same 
spatial resolution as Kuo et al. (2008) for their active 
galaxy sample. Because the control sample is a factor 
of two closer in distance, their HI intensity is a factor of 
4 higher for the same gas mass. Everything else being 
equal, this would require a factor of 16 shorter integra- 
tion time for our control sample, a considerable saving 
given that the in tegration time for each active galaxy in 



Kuo et al. (2008) was ^2 hrs. In practise, to achieve the 
same spatial resolution of ^20 kpc (i.e., angular resolu- 
tion of ~120'') as in the active galaxy sample (which was 
imaged at an angular resolution of ~60''), we had to dis- 
card baselines that were more than one-third the length 
of the longest baselines in the active galaxy sample (i.e., 
using baseline up to only 1.5 kA). Because there are 
more shorter than longer baselines, the actual integra- 
tion time per object in our control sample was therefore 
~12 mins (still a considerable saving in time over 2 hrs) 
and the observing time per object ~20 mins (compared 
with ~2.5 hrs for the active galaxy sample). In this way, 
we achieved a rms fiux density of ~1.8 mJy beam~^ in 
one channel, corresponding to a 5cr detection threshold in 
HI gas column density of ~3 x 10^^ cm~^ per synthesized 
beam of diameter ~120'' (20 kpc) in a given channel with 
a velocity width ~25.2 km s"-*^. This is the same detec- 



tion threshold as reached by |Kuo et al.| ([2008| in their 
active galaxy sam ple. 



Like Kuo et al. (2008), we performed all the data re- 
duct ionTE'TIie'stanaixa manner using the NRAO AIPS 
package. Unlike for the active galaxy sample, we found 
no radio frequency interference (RFI) in the data for our 
control sample, presumably because their redshifted HI 
lines lie closer to the protected band at 21 cm. Fol- 
lowing bandpass calibration, as well as time-dependent 
amplitude and phase calibration, we subtracted the con- 
tinuum emission from the visibility data by interpolating 
between line-free channels (determined through trial and 
error) on either side of the HI line. For NGC 4256, where 



there is a relatively strong continuum source in the field, 
we first subtracted the continuum emission of that source 
in the visibility plane. Finally, we made maps of each 
channel, and corrected for the primary beam response 
of the antennas. We also combined the channel maps 
(before primary beam correction) to make maps in total 
intensity (zeroth moment) and intensity-weighted mean 
velocity (first moment), where we tried several different 
combinations of spectral smoothing and minimum signal 
cutoff level in each channel so as to suppress the noise 
and bring out faint diffuse features. 

4. RESULTS 

We detected twenty-one of the twenty-seven galaxies 
in our control sample. Two of the six galaxies not de- 
tected, NGC 4128 and NGC 4543, have (to the best of 
our knowledge) not previously been observed in HI. The 
reported upper limits in the integrated HI intensities for 
NGC 4578 (Hutchmeier et al. 1998) and NGC 5326 
(Theureau et al. 1998) are below that attained in our 
observations. NGC 3468 has a reported integrated HI in- 
te nsity of 1.40 Jy km s ~^ (with no reported uncertainty) 
in 



Huchtmeier 



(1989), and should therefore have been 

detected in our observatio n. On checking the o rigina l 
source of the detection in Giovanelli & Haynes' ("1981"), 
we found that the actual galaxy det ected was UGC 346 8 
and not NGC 3468 as reported in Huchtmeier (1989). 
NGC 5355 has a repo rted integrat e d HI intensity of 



14.6 ± 1.8 Jy km s'^ in [Huchtmeier | ( |1989[ ), and should 
therefore have been detected in our observation. On 
checking the original source of the detection, Richter and 
Huchtmeier (1991) mentioned that the 8.^8 beam of the 
Effelsburg telescope used for this observation would in- 
clude also NGC 5353 and NGC 5354. We note that this 
beam also includes NGC 5350, which is one of the galax- 
ies in our control sample. We detected NGC 5350 at a 
much higher integrated HI intensity than that reported 
by Richter and Huchtmeier (1991) for NGC 5355, sug- 
gesting that the actual HI detection was of NGC 5350. 

The twenty-one galaxies that we detected in HI com- 
prise our ensemble control sample from which we draw 
statistical results. For comparison, there a re twenty- 
three a ctive galaxies in the ensemble sample of [Kuo et al. 
(2008), of which eighteen are classified as Seyfert galax- 
ies] The HI maps of the ensemble control sample are 
shown in Figures 2-25, and their HI properties summa- 
rized in Table 3. We have overlaid each integrated HI 
intensity (zeroth moment) map on an optical (either the 
red or blue filters) image from the 2nd Digitized Sky Sur- 
vey (DSS2). This database provides a relatively uniform 
set of optical images to search for any disturbances in 
optical starlight. 

We also detected in HI a number of galaxies in the 
same field of view as our target galaxies. Those galax- 
ies lying within the primary beam (i.e., within a radius 
of ~ 15.^8 about phase center) are listed in Table 4, and 
are referred to by their names from optical catalogs. We 
also list their optical redshifts where available, as well as 
their corrected apparent 5-band magnitudes and Hubble 
types as listed in the Hyperleda database. The measured 
HI properties of these neighboring galaxies are summa- 
rized in Table 5. Apart from neighboring galaxies de- 
tected in HI, we also searched the Hyperleda database in 
a radius of ~15' about each galaxy in our control sample 
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to obtain a measure of the richness of their fields. We 
emphasize that this is by no means an exhaustive search. 

4.1. Individual Galaxies 

We present here the results for each galaxy, highl ight 
ing features of most relevance to this work. As in Kuo 



et al. (2008), we separate the individual galaxies in our 
ensemble control sample into the following three groups. 

1. Those that clearly show HI tidal features tracing 
interactions with neighboring galaxies. These fea- 
tures are in the form of tidal bridges that connect 
the two interacting galaxies, or an extension from 
one or both galaxies in the direction of the other 
(i.e., an incomplete tidal bridge). In addition, fea- 
tures in the form of tidal tails comprising a pro- 
trusion or curved extension from one galaxy on the 
side away from the other interacting galaxy may 
be seen. Only one of the twenty-one galaxies in 
our ensemble control sample of inactive galaxies 
fall into this group. By comparison, thirteen of 
the twenty-three galaxies in the ensemble sample 
of active galaxies, including twelve of the eighteen 
classified as Seyferts, fall into this group. 

2. Those that clearly show both spatial and kinematic 
disturbances in HI, but which cannot be directly 
linked to interactions with neighboring galaxies if 
any. These disturbances are in general less promi- 
nent than the tidal features observed in the first 
group. Only two of the twenty-one galaxies in our 
ensemble control sample of inactive galaxies fall 
into this group. Four of the twenty-three galaxies 
in the ensemble sample of active galaxies, including 
three of the eighteen classified as Seyferts, fall into 
this group. 

3. Those that show marginal or no detectable HI dis- 
turbances. Eighteen of the twenty-one galaxies in 
our ensemble control sample of inactive galaxies fall 
into this group. Only one exhibits weak HI dis- 
turbances; the remaining seventeen exhibit no dis- 
cernible HI disturbances. By comparison, only six 
of the twenty-three galaxies in the ensemble sam- 
ple of active galaxies, including only three (two of 
which exhibit weak HI disturbances) of the eigh- 
teen classified as Seyferts, fall into this group. 

Because the vast majority of the galaxies in our control 
sample exhibit no detectable HI disturbances at an angu- 
lar resolution of ~120'' (~20 kpc), for these undisturbed 
galaxies we show only HI maps at the full angular reso- 
lution attained in our observations of ~60'' (~10 kpc). 
The maps at ~120'' angular resolution contain no ex- 
tra information. For the disturbed galaxies, moment 
maps are shown at angular resolutions of both ~120'' 
and ~60'' for the reader to compare, and HI channel 
maps at an angular resolution of ~120'' only. The HI 
channel maps of the undisturbed galaxies are available in 
the electronic version of this paper (Figs. 28-44). (For as- 
troph readers, the channel maps are available at the link 
|http://www. astro. virginia.edu~ck2v/maps) At an angu- 
lar resolution of ~60'', the corresponding 5cr detection 
threshold in HI gas column density is ~7 x 10^^ cm~^ 



per synthesized beam in a given channel with a velocity 
width ~25.2 km s~^, which is higher than the Seyfert 
sample. Because the Seyfert sample is twice further away 
from us comparing with the control sample, the detec- 
tion limit in terms of the HI mass of the control sample 
(~ 6 xlO^ Mq) is lower than of the seyfert sample (^ 7 
xlO^Mo). 

4.1.1. Group I: HI tidal features tracing interactions 
with neighboring galaxies 

At a spatial resolution of ~20 kpc only one of our 
ensemble sample of inactive galaxies exhibit tentative 
HI disturbances, which only in the map at a spatial 
resolution of 10 kpc can be clearly seen to trace to tidal 
interactions with neighboring galaxies. Nevertheless, 
we generously place this sole example of its kind in 
Group 1. When categorizi ng the ensemble sample of 
active galaxies in Kuo et al.| ( |2008[ ), we took a more 
conservative approach. 

NGC 2543 

NGC 2543 is a barred Sb galaxy (all Hubble types from 
the Hyperleda database, just as for the active galaxies in 
our companion paper). It appears somewhat asymmet- 
ric in the optical, although not obviously disturbed. A 
relatively dim galaxy at a comparable optical redshift, 
PGC 80408, lies -3f4 (-36 kpc) south of NGC 2543. 

Our HI image at a spatial resolution of —20 kpc (Fig. 2) 
reveals that NGC 2543 shares a common envelope with 
PGC 80408. An examination of the channel maps at this 
spatial resolution (Fig. 3) reveals tentative evidence for a 
tidal bridge (which, if real, is at best marginally spatially 
resolved) connecting the two galaxies. Our HI image at a 
spatial resolution of —10 kpc (also shown in Fig. 2), and 
also the channel maps at this angular resolution (also 
Fig. 3), clearly reveals that the connection between the 
two galaxies is indeed a tidal bridge. In addition, the 
HI disk of NGC 2543 is clearly disturbed on its south- 
western side, which on this side extends much further 
than on the north-eastern side and curls in a direction 
towards its interacting neighbour PGC 80408. 

4.1.2. Group II: HI disturbances likely produced by 
interactions 

Two of the galaxies in our ensemble control sample 
exhibit spatial and kinematic disturbances in HI, but 
unlike the situation in Group I the observed perturba- 
tions cannot be immediately linked to interactions with 
neighboring galaxies. 

NGC 2551 

NGC 2551 is a Sa galaxy. A dimmer galaxy at a compa- 
rable optical redshift, UGC 4390, lies -14;4 (-146 kpc) 
to the north-east of NGC 2551. Neither galaxies appear 
to be visibly disturbed in the optical. 

We detected both NGC 2551 and UGC 4390 in HI, as 
well as the relatively faint galaxy PGC 2755603, which 
has no previously reported optical redshift (Figs. 4 and 
5). At a spatial resolution of —20 kpc (Fig. 5), the HI 
disk of NGC 2551 appears to be disturbed, with this 
disturbance especially prominent in the first moment 
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map (kinematics). The observed disturbance in the HI 
disk is even more prominent at an angular resolution 
of ~10 kpc (Fig. 4), but not easily apparent in the 
channel maps (Fig. 6). By comparison, UGC 4390 
does not appear to be disturbed in HI at either spatial 
resolutions, although it is possible that the HI kinematic 
axis is not aligned with the optical major axis (the 
optical major axis is poorly defined in the outer regions 
of this galaxy). PGC 2755603 was not spatially resolved 
in HI. It is not clear whether the observed disturbance 
in NGC 2551 was caused by tidal interactions with (one 
of) these or other galaxies. 

NGC 3094 

NGC 3094 is a barred Sa galaxy that does not appear 
to be obviously disturbed in the optical. In HI (Fig. 7), 
we detected NGC 3094 and also PGC 2806972, a rela- 
tively faint galaxy that liest ~7f5 (~78 kpc) to the east 
of NGC 3094. PGC 2806972 has no previously reported 
optical redshift. 

NGC 3094 exhibits a HI protrusion on the south- 
eastern side that is kinematically discontinuous from 
its HI disk. A close examination of the channel maps 
(Fig. 8) reveals that this protrusion in physically linked 
to NGC 3904. These channel maps also reveal prominent 
north-south extensions over a narrow velocity range not 
seen in the moment maps; these extensions are orthog- 
onal to the HI kinematic axis of the galaxy. There are 
no cataloged galaxies at the location of this protrusion, 
nor any optical counterpart visible in the DSS2 image. 
It is possible that this protrusion is an incomplete tidal 
bridge caused by interactions with PGC 2806972. 

4.1.3. Group III: Weak or no HI disturbances 

Eighteen of the twenty-one galaxies fall into this group, 
accounting for the vast majority of the galaxies in our 
ensemble control sample. Of these eighteen, seventeen 
do not show any evidence whatsover for disturbances 
in HI. Only NGC 3976, which as we explain below 
has recently been classified as an active galaxy, shows 
evidence for weak disturbances. 

NGC 3976 

NGC 3976 is a barred Sb galaxy. It is clearly dis- 
turbed in the optical (Fig. 9), exhibiting asymmetric ex- 
tensions on the north-eastern and south-western sides of 
the disk. Alth ough NGC 3976 is not classifi ed as an ac- 



(1998), the source 



tive galaxy by Veron-Cetty & Veron 
used to select the active galaxy sample, it is classified as 
a Seyfert 2 with Mb = —21.1 in the latest version of this 
catalog (Veron-Cetty fc Veron,, 2006) . Thus, if NGC 3976 
was placed at the same range of redshifts as the active 
galaxy sample, it would have been selected as part of 
that sample. For internal consistency, however, we re- 
tain NGC 3976 in our control sample. If included in the 
active galaxy sample, this would make the conclusions 
reached in our study only stronger. A relatively faint 
galaxy at a comparable optical redshift, PGC 37490, lies 
-4f5 to the south of NGC 3976. 

We detected NGC 3976, but not PGC 37490. At a 
spatial resolution of ~20 kpc, the HI disk of NGC 3976 
appears to be asymmetric (Fig. 9). The HI gas extends 



further out on the north-eastern side compared to the 
south-western side, giving the disk a lopsided appearence 
in HI. In addition, the HI kinematic axis appears to lie at 
a different position angle compared with the major axis 
of the optical disk. At a spatial resolution of ~10 kpc 
(also shown in Fig. 9), this asymmetric appearence 
appears to be caused in part (but perhaps not wholly) 
by the larger HI extension of the north-eastern compared 
with the south-western spiral arm. 

IC 2461 

IC 2461 is an edge-on Sb galaxy that does not appear 
to be visibly disturbed in the optical. There are no cat- 
aloged galaxies at a comparable optical redshift within 
15' of the galaxy. 

At a spatial resolution of ~20 kpc, the HI disk of 
IC 2461 appears to be symmetric both spatially and 
kinematically with no evidence for perturbations. At a 
spatial resolution of ~10 kpc (Fig. 10), the HI major 
and kinematic axes appears to exhibit an aticlockwise 
twist on the north-western side of the disk, perhaps 
indicative of a warp. 

NGC 3835 

NGC 3835 is a Sab galaxy that is not visibly disturbed 
in the optical. There are no cataloged galaxies at a com- 
parable optical redshift within 15' of the galaxy. 

At a spatial resolution of ~20 kpc, the HI disk of 
NGC 3835 appears to be symmetric both spatially and 
kinematically with no evidence for perturbations. The 
same also is true at a spatial resolution of ~10 kpc 
(Fig. 11), where the central inward contraction of 
the intensity contours in the zeroth moment maps is 
probably caused by a depression in HI intensity at the 
center of the galaxy (as is often observed in nearly disk 
galaxies). 

NGC 4067 

NGC 4067 is a Sb galaxy that is not visibly disturbed 
in the optical. There are no cataloged galaxies at a com- 
parable optical redshift within 15' of the galaxy. 

At a spatial resolution of ~20 kpc, the HI disk of 
NGC 4067 appears to be symmetric both spatially 
and kinematically with no evidence for perturbations. 
At a spatial resolution of ~10 kpc (Fig. 12), the HI 
disk appears to be more prominent or extended on the 
south-western side compared with the north-eastern 
side. Other than this small asymmetry, the HI disk does 
not exhibit any spatial nor kinematic disturbances. Once 
again, the central inward contraction of the intensity 
contours in the zeroth moment maps is probably caused 
by a depression in HI intensity at the center of the galaxy. 

NGC 4256 

NGC 4256 is an edge-on Sb galaxy that is not visi- 
bly disturbed in the optical. There are five relatively 
faint cataloged galaxies at a comparable optical redshift 
within 15' of NGC 4256. Their separations range from 
0;4 (-4 kpc) to 14' (-153 kpc) from NGC 4256. None 
of these neighboring galaxies were detected in HI. 
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The HI gas mass of NGC 4256 is the lowest among 
all of the galaxies in our ensemble control sample, and 
hence the corresponding signal-to-noise ratio of the HI 
map for this galaxy poorest. At a spatial resolution 
of -20 kpc, the HI disk of NGC 4256 appears to be 
symmetric both spatially and kinematically with no 
evidence for perturbations. The same also is true at a 
spatial resolution of —10 kpc (Fig. 13). 

NGC 4275 

NGC 4275 is a Sb galaxy that has the smallest incli- 
nation among all of the galaxies in our ensemble control 
sample. There are no cataloged galaxies at a comparable 
optical redshift within 15' of the galaxy. 

At a spatial resolution of —20 kpc, the HI disk of 
NGC 4275 appears to be symmetric both spatially and 
kinematically with no evidence for perturbations. The 
same also is true at a spatial resolution of —10 kpc 
(Fig. 14). 

NGC 4351 

NGC 4351 is a barred Sab type galaxy that appears 
to be disturbed in the sense that it is lopsided in the 
optical, being more extended on the western than eastern 
side of the galaxy. There are no cataloged galaxies at a 
comparable optical redshift within 15' of the galaxy. 

At a spatial resolution of —20 kpc, the HI disk of 
NGC 4351 appears to be symmetric both spatially and 
kinematically with no evidence for perturbations. At a 
spatial resolution of —10 kpc (Fig. 15), the HI disk is 
asymmetric in that it is more prominent or extended on 
the western side, just like in the optical. 

NGC 4384 

NGC 4384 is a Sa galaxy that is not visibly dis- 
turbed in the optical. There is a relatively faint galaxy, 
PGC 2478952, with a comparable optical redshift at 
-14.'9 (-162 kpc) from NGC 4384. That galaxy was 
not detected in HI. 

At a spatial resolution of —20 kpc, the HI disk of 
NGC 4384 appears to be symmetric both spatially and 
kinematically with no evidence for perturbations. At 
a spatial resolution of —10 kpc (Fig. 16), the HI disk 
may be somewhat asymmetric in that the eastern side 
appears to be more prominent or extended. 

NGC 4470 

NGC 4470 is a Sa type galaxy that is not visibly dis- 
turbed in the optical. The galaxy that is closest in optical 
redshift and which lies within 15' of NGC 4470 differs in 
systemic velocity by —900 km s~^, and lies outside the 
velocity coverage of our HI observation. 

At a spatial resolution of —20 kpc, the HI disk of 
NGC 4470 appears to be symmetric both spatially and 
kinematically with no evidence for perturbations. At 
a spatial resolution of —10 kpc, the HI disk may show 
weak asymmetries (Fig. 17). 

NGC 4513 



NGC 4513 is a SO galaxy that is not visibly disturbed 
in the optical. There are no cataloged galaxies at a com- 
parable optical redshift within 15' of the galaxy. 

At a spatial resolution of —20 kpc, the HI disk of 
NGC 4513 appears to be symmetric both spatially and 
kinematically with no evidence for perturbations. The 
same also is true at a spatial resolution of —10 kpc 
(Fig. 18). This galaxy has the second lowest HI gas mass 
among our ensemble control sample, with the resulting 
low signal-to-noise ratio making the presumed central 
depression in HI gas particularly prominent. 

NGC 4567 and NGC 4568 

Both NGC 4567 and NGC 4568 are Sbc galaxies that 
form an optically-overlapping pair with comparable op- 
tical redshifts. In this way, NGC 4567/NGC 4568 resem- 
ble UGC 3995/UGC 3995A in the active galaxy sample. 
Both NGC 4567 and NGC 4568 are included in our con- 
trol sample. There are no other cataloged galaxies at a 
comparable optical redshift within 15' of the galaxy. 

Despite their apparent proximity and similar redshifts, 
as in the case of UGC 3995/UGC 3995 A we detect no 
HI disturbances at a spatial resolution of —20 kpc, nor 
at a spatial resolution of —10 kpc (Fig. 19), in either 
NGC 4567 or NGC 4568. Indeed, in a much deeper HI 
observation at an angular resolution of —20", |Iono et al. 
( ,2005) detected no extended tidal features in this pair 
of galaxies. 

NGC 4591 

NGC 4591 is a Sb galaxy that is not visibly disturbed 
in the optical. There are no cataloged galaxies at a com- 
parable optical redshift within 15' of the galaxy. 

At a spatial resolution of —20 kpc, the HI disk of 
NGC 4591 appears to be symmetric both spatially and 
kinematically with no evidence for perturbations. The 
same also is true at a spatial resolution of —10 kpc 
(Fig. 20). 

NGC 4814 

NGC 4814 is a Sb galaxy that is not visibly dis- 
turbed in the optical. There is a relatively faint galaxy, 
SDSSJ 125426.85+582348.9, with a comparable optical 
redshift at -7.'9 (-86 kpc) from NGC 4384. That galaxy 
was not detected in HI. 

At a spatial resolution of —20 kpc, the HI disk of 
NGC 4814 appears to be symmetric both spatially and 
kinematically with no evidence for perturbations. The 
same also is true at a spatial resolution of —10 kpc 
(Fig. 21). 

NGC 4964 

NGC 4964 is a SO-a galaxy that is not visibly disturbed 
in the optical. There are no cataloged galaxies at a com- 
parable optical redshift within 15' of the galaxy. 

At a spatial resolution of —20 kpc, the HI disk of 
NGC 4964 appears to be symmetric both spatially and 
kinematically with no evidence for perturbations. The 
same also is true at a spatial resolution of —10 kpc 
(Fig. 22). 
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NGC 5289 

NGC 5289 is a nearly edge-on barred Sa galaxy that 
is not visibly disturbed in the optical. There are six 
mostly relatively faint cataloged galaxies at a compara- 
ble optical redshift within 15' of NGC 5289. Their sep- 
arations range from 5f3 (~58 kpc) to 12'7 (~138 kpc) 
from NGC 5289. Only of these neighboring galaxies, 
NGC 5290, lying furthest away from and brighter in the 
optical than NGC 5289, was detected in HI. 

At a spatial resolution of ~20 kpc, the HI disk of 
NGC 5289 appears to be symmetric both spatially and 
kinematically with no evidence for perturbations. The 
same also is true at a spatial resolution of ~10 kpc 
(Fig. 23). 

NGC 5350 

NGC 5350 is a barred Sbc galaxy that appears to be 
somewhat asymmetric although not obviously disturbed 
in the optical. There are seven cataloged galaxies at 
a comparable optical redshift within 15' of NGC 5350, 
two of which are comparably bright and the remainder 
relatively faint. None of these neighboring galaxies were 
detected in HI. 

At a spatial resolution of ~20 kpc, the HI disk of 
NGC 5350 appears to be symmetric both spatially and 
kinematically with no evidence for perturbations. The 
same also is true at a spatial resolution of ~10 kpc 
(Fig. 24). 

NGC 5375 

NGC 5375 is a barred Sab galaxy that is not visibly dis- 
turbed in the optical. There is a relatively faint galaxy, 
PGC 49623, with a comparable optical redshift at ~4.'6 
(~47 kpc) from NGC 5375. That galaxy was not de- 
tected in HI. Instead, we detected the galaxy MAPS- 
NGP 0-325-1536477 (as cataloged in the NASA Extra- 
galactic Databse), which has no previously reported op- 
tical redshift, which is 9'A (~131 kpc) to the south-east 
of NGC 5375. 

At a spatial resolution of ~20 kpc, the HI disk of 
NGC 5375 appears to be symmetric both spatially and 
kinematically with no evidence for perturbations. The 
same also is true at a spatial resolution of ~10 kpc 
(Fig. 25). 

4.2. Other Galaxies in Target Fields 

Apart from the optically-overlapping pair 
NGC 4567/NGC 4568, we detected only six other 
galaxies in HI within the primary beam of our target 
fields. None of these galaxies exhibit HI disturbances at 
either 20 kpc or 10 kpc spatial resolution. The HI maps 
of these galaxies are shown together with the maps for 
the control sample galaxies. 

4.3. Ensemble Statistics 

We now make statistical inferences from the results 
based on our ensemble control sample of twenty-one 
galaxies that we mapped in HI. Only four of the twenty- 
one (~19%) galaxies in the control sample exhibit de- 
tectable disturbances at a spatial resolution of ~20 kpc. 



By contrast, seventeen of the twenty-one (~81%) galax- 
ies exhibit no HI disturbances whatsoever at this spatial 
resolution. One of these disturbed galaxies (NGC 3976) 
has since been found to be a Seyfert 2 galaxy, and so only 
three of the twenty (15%) inactive galaxies are actually 
disturbed in HI, whereas seventeen of twenty (~85%) 
show no HI disturbances. The HI disturbances be di- 
rectly traced to tidal interactions with neighboring galax- 
ies in only one case. In the remaining two cases with HI 
disturbances, the observed disturbances are likely to be 
produced by tidal interactions with neighboring galax- 
ies also detected in HI. By comparison, seventeen of the 
eighte en (^94%) Seyfer t galaxies in the ensemble sam- 



ple of Kuo et al. (2008) exhibit HI disturbances, twelve 
(^67%) of which are mteracting with neighboring galax- 
ies (Group I), another three (~17%) probably interacting 
with neighboring galaxies (Group II), and the remaining 
two (~12%) weakly disturbed (Group HI). 

Only two of the twenty-one (~ 10%) galaxies in the 
control sample are visibly disturbed in optical DSS2 im- 
ages. Leaving aside NGC 3976 since found to be a 
Seyfert 2 galaxy, only one of twenty (~ 5%) inactive 
galaxies are optically disturbed. By comparison, six of 
the twenty-three (^26%) galaxies in the ensemble s am- 
ple of active galaxies studied by Kuo et al. (2008) are 
visibly disturbed in optical DSS2 images, mcluding five 
of the eighteen galaxies (~28%) classified as Seyferts. 
As pointed out in § 1, the selection criteria used may 
select against (strongly) optically disturbed galaxies in 
the control sample. Nonetheless, the relatively small dif- 
ference in the frequency of optical disturbances between 
the Seyfert and control samples (~28% vs ~5%) pales 
besides the dramatic difference in the frequency of HI dis- 
turbances between these two samples (~94% vs ~15%). 

In case of any severe selection bias against optically- 
disturbed galaxies in our control sample, here we com- 
pile statistics for just the optically-undisturbed galaxies 
in both the Seyfert and control sample. Of the thir- 
teen such Seyfert galaxies, twelve (~92%) exhibit HI 
disturbances, nine (~69%) of which are interacting with 
neighboring galaxies (Group I), one (~7%) possibly in- 
teracting with a neighboring galaxy (Group II), and two 
(~15%) weakly disturbed (Group III). By comparison, 
only three of the twenty (~15%) optically-undisturbed 
inactive galaxies are disturbed in HI. 

5. INTERPRETATION AND DISCUSSION 

5.1. Relative HI Gas Masses 

We first examine whether the dramatic difference in 
HI disturbances between the active and control samples 
might be caused by differences in their HI gas masses. 
This could happen if, for example, the same fractional 
HI gas mass is displaced into tidal features of equal spa- 
tial dimensions in a given interaction. In such a case, 
interacting galaxies with lower HI gas masses will ex- 
hibit dimmer tidal features, which may fall below the 
detection threshold in a given observation. 

In Figure 26, we plot the HI gas masses of the active 
and control samples. As can be seen, the active sample is 
distributed towards somewhat higher HI gas masses than 
the control sample. A K-S test gives a statistic D = 0.30 
(converging to for identical distributions) and prob- 
ability Pr = 0.23 (1 for identical distributions). More 
importantly, most of the galaxies in our control sample 
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span the same range of HI gas masses as the active sam- 
ple for which we detect HI disturbances. We therefore 
believe that the dramatic difference in HI disturbances 
between the two samples cannot be caused by the rel- 
atively small difference in their overall distributions of 
HI gas masses. The possibility that tidal features in our 
control sample somehow conspire to be less spatially ex- 
tended than those in the active sample is ruled out by our 
HI images of the control sample at a spatial resolution of 
~10 kpc. 

5.2. Incidence of Neighboring Galaxies 

In Figure 27, we show the cumu lative fraction of 
Seyfert galaxies (solid line) studied bylKuo et al.j( [2QQ8 ) 
with (candidate) interacting neighboring galaxies (i.e., in 
Groups I, II, or HI), plotted as a function of their pro- 
jected separations (for the few with multiple interact- 
ing neighboring galaxies, their nea rest such neighbor). 
As explained in |Kuo et al. (2008), over the range of 
projected separations plotted (up to 90 kpc), there are 
only two neighboring galaxies not identified to be (possi- 
bly) interacting with their respective (the same) Seyfert 
galaxy, the third of which is identified as the interact- 
ing neighbor. Thus, including all possible neighboring 
galaxies around Seyfert galaxies do not change this fig- 
ure, at least within the range of projected separations 
plotted. In the same figure, we also plot the cumulative 
fraction of galaxies in our control sample (dashed line) 
with neighboring galaxies within the same range of pro- 
jected separations. (Keep in mind that our observations 
span the same velocity interval for both samples.) Our 
detection threshold is relatively uniform (within ~20%) 
and essentially the same for both the Seyfert and control 
samples at projected separations up to ~90 kpc. As can 
be seen, the frequency of Seyfert with neighboring galax- 
ies is clearly much higher than that of the control sample 
over the range of projected separations plotted. This dif- 
ference simply reflects the prevalence of tidal interactions 
in the Seyfert but not control sample. 

As mentioned in §1, some optical studies find an excess 
of Seyfert with (projected) neighboring galaxies com- 
pared with inactive galaxies, whereas others do not. A 
discussion of the different methodologies used in differ- 
ent studies, as well as reasons for their conflicting re- 
sults, is beyond the scope of this paper. Any compar- 
ison made between previous optical and our HI studies 
should keep in mind that our HI imaging study is not 
designed to address differences between the fraction of 
Seyfert and inactive galaxies with neighboring galaxies. 
(Rather than looking for any such difference, we directly 
determine what fraction of Seyfert and inactive galax- 
ies show HI disturbances, and where possible determine 
whether these disturbances are caused by tidal int erac- 
tions with n eighboring galaxies.) Both the study of Kuo 
et al. (2008) and that reported here can only detect rela- 
tively gas-rich galaxies, and so a HI census of neighboring 
galaxies is likely to be less complete than an optical cen- 
sus. 

Nevertheless, to see what can be learnt, we compare 
th e results of our HI stu dies with t he optical studies 



of [Dultzin-Hac yan et al.| (1999) and |Koulouridis et al. 
( 2006|) , who measured the fraction of galaxies in both 
their h^eyfert and control samples with (projected) neigh- 
boring galaxies as a function of their projected separa- 



tions. Both these studies employed the same Seyfert and 
control sa mples; the Seyfert sam ple was taken from the 
catalog of Lipovetsky et al. ( 1988"), which contained most 
if not all of the then known Seyfert galaxie s, and the 



ma tched control sample from the Cf A survey ( Huchra et 

[iI ||T983l). 

Dultzin-Hacyan et al. (1999) searched for projected 
neighboring galaxies withm a radius of ~140 kpc us- 
ing the (first) Digitized Sky Survey (DSS). They found 
that in both their Seyfert and control samples, the frac- 
tion having projected neighboring galaxies with diam- 
eters 4-10 kpc increases as the search radius widens, 
reaching between 80% -100% within a search radiu s of 

At 



1999) 



■^100 kpc (see Fig. 1 of Dultzin-Hacyan et al.[ 

any given search radius, the fraction of Seyfert 1 and in- 
active galaxies possessing projected neighboring galaxies 
is comparable. On the other hand, there is a relatively 
large difference between the fraction of Seyfert 2 and in- 
active galaxies possessing projected neighboring galaxies 
at the smallest search radius of 10 kpc, with this differ- 
ence diminishing as the search radius widens until dis- 
appearing at ~100 kpc. Thus, the difference, if any, be- 
tween the fraction of Seyfert 2 and matched samples of 
inactive field galaxies possessing projected neighboring 
galaxies depends at least in part on the search radius 
used. We note that this may explain the seemingly con- 
tradictry results between some optical studies; i.e., those 
using relatively small search radii finding a higher frac- 
tion of Seyfert with projected neighboring galaxies than 
in matched samples of inactive galaxies, whereas those 
using relatively large search radii finding no appreciable 
difference. 



Our results agree with those of [Dultzin-Hacyan et al. 
(1999) in two important r espects. First, nearly all the 
Seyfert galaxies studied by 'Dultzin-Hacyan et al. (1999) 
have projected neighboring galaxies within ~ 100 Kpc. 
Second, we both see an excess of Seyfert (our ensemble 
sample contains a much larger fraction of Seyfert 2 than 
Seyfert 1) having (projected) neighboring galaxies at pro- 
jected separations smaller than ~90 kpc compared with 
the control sample. Where we differ is in our seeing this 
excess becoming increasingl y larger with projected sepa- 
ration up to 90 kpc, whereas Dultzin-Hacyan et al. ( 1999 ) 
find this excess to become increasingly smaller witn pro- 
jected separation before disappearing at ~100 kpc. Keep 
in mind, however, that our HI imaging observations 
cleanly pick out genuine neighboring but only gas-rich 



galaxies, whereas those of Dultzin-Hacyan et al.| (1999) 
may include spurious obje cts. 

'Koulouridis et al.l( 2006 ) repeated the study of'Dult ziiT] 
Hacyan et al. (1999| but now armed with redshifts from 
the second CfA as well as Southern Sky Redshift Sur- 
veys. Confi r ming the trends s e en by ,Dultzin-Hacyan| 
et al.| (|1999D , [Koulouridis et al.| ( [2006[ ) find a relatively 



large difference between the fraction of Seyfert 2 (but, 
again, not Seyfert 1) with neighboring galaxies within a 
velocity difference of 200 km s~^ compared with inactive 
galaxies at the smallest measured projected separation of 
7 kpc, with this difference diminishing as the search ra- 
dius widens until disappearing at roughly 70 kpc (scaled 
to Ho = 67 km s~^ Mpc~^). With a maximum search 
radius of ~70 kpc and a limiting magnitude of ^ 15.5 
(1 magnitude brighter than the LMC at the distance of 
our sample), this study would have picked up only a small 
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fraction of th e interacting neighboring galaxies that Kuo 
!008|) detected around t heir Seyfert sample. In- 
oulouridis et al. (2006) find that only 



deed, 

Seyfert 1 and 



. 10% of 

30% of Seyfert 2 galaxies have neigh- 



boring galaxies within the range of parameters searched. 
Going to a lower limiting magnitude of ~ 18.5 (cor 
responding to the SMC at th e distance of our sample 
for a subset of their sample, [Koulouridis et alT ( |2006 
find that about twice as many still of both Seyfert 
and Seyfert 2 galaxies have neighboring galaxies within 
a projected separation of ~50 kpc. This fraction is com- 
parable to that we found in our Seyfert sample (~40% 
at a projected se 



paration of ~50 kpc). Unfortunately, 
20061) did not search for neighboring 



Koulouridis et al 
galaxies down to the the same low limiting magnitudes 
around their control samples, and hence a more detailed 
comparison cannot be made. 



5.3. Implicating Tidal Interactions in Triggering AGNs 

The contrast between the prevalence of HI distur- 
bances in the active and control samples is dramatic. In 
the active galaxy sample, we can only find a few galaxies 
that are not disturbed spatially and usually also kinemat- 
ically on galactic- wide (> 20 kpc) scales. In the control 
sample, we find the opposite result, with most exhibiting 
no detectable spatial or kinematic disturbances on both 
the same as well as smaller (> 10 kpc) scales. 

Specifically, of the eighteen galaxies in the ensem- 
ble sample of twenty-three active galaxies classified as 
Seyferts, seventeen (~94%) exhibit HI disturbances. In 
at least ~67% and possibly as high as ~94% of cases, the 
observed HI disturbances can be traced to tidal interac- 
tions with neighboring galaxies detected also in HI. By 
contrast, only four of the twenty-one (~19%) galaxies in 
our control sample exhibit HI disturbances. Removing 
the one disturbed galaxy since found to be a Seyfert 2 
galaxy, only three of the twenty (15%) are actually dis- 
turbed in HI. Only one of the eighteen (~6%) Seyfert 
galaxies in the ensemble active galaxy sample exhibit no 
detectable HI disturbances whatsoever. By contrast, sev- 
enteen of the twenty (~85%) galaxies in our control sam- 
ple exhibit no detectable HI disturbances. These results 
directly implicate tidal interactions in initiating events 
that lead to optically-luminous Seyfert activity in a large 
fraction of local disk galaxies. 

6. SUMMARY AND CONCLUSIONS 

The central purpose of this paper is to determine 
whether the high incidence of tidal interactions observed 
in HI gas for a sample of Seyfert galaxie s reported in 
the companion paper by Kuo et al. (2008) is related to 
their AGN activity. Our strategy was to image at the 
same spatial resolution (~20 kpc) and sensitivity in HI 
gas a comparably large number of inactive galaxies that 
were closely matched in Hubble type and to the degree 
possible optical luminosity, as well as range in size and 
inclination, to the active galaxies. We detected twenty- 
one of the twenty-seven galaxies in our control sample, 
imaged at the same spatial resolution and HI column den- 
sity threshold as the active galaxy sample. These twenty- 
one galaxies comprised our ensemble control sample from 
which we drew the following statistical results. 

1. Only four of the twenty-one galaxies (~19%) ex- 
hibit spatial and usually also kinematic distur- 



bances on galactic (> 20 kpc) scales. One of these 
disturbed galaxies has since been found to be a 
Seyfert 2 galaxy, and so only three of the twenty 
(15%) inactive galaxies in our ensemble control 
sample are actually disturbed in HI. 



2. Seventeen of the twenty-one galaxies (~81%) show 
no HI disturbances whatsover on galactic (> 
20 kpc) scales. Excluding again the one disturbed 
galaxy since found to be a Seyfert 2 galaxy, seven- 
teen of the twenty (85%) inactive galaxies in our 
ensemble control sample are not disturbed in HI. 

By contrast, of the eighteen galaxies in the ensem- 
ble sample of twenty-three active galaxies classified as 
Seyferts, seventeen (~94%) exhibit HI disturbances. In 
at least ~67% and possibly as high as ~94% of cases, 
the observed HI disturbances can be traced to tidal in- 
teractions with neighboring galaxies detected also in HI. 
The dramatic contrast in the incidence of HI disturbances 
between the active and inactive galaxy samples strongly 
implicates tidal interactions in initiating events that lead 
to luminous Seyfert activity in a large fraction of local 
disk galaxies. 
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Fig. 1. — Number distribution of the active galaxy sample selected by \Kuo et ah] [200 87 (unfilled histogram) and the inactive control 
sample selected here (hatched histogram) in morphologial type (upper panel j, corrected absolute B-band magnitude (second panel from top, 
M^), optical size (second panel from bottom), and inclination (bottom panel). The morphological type, M^, optical size and inclination 
angle are retrieved from the HYPERLEDA database. Note that there are 27 active samples included in the plot. The of MS 
04595+0327 is not available in the HYPERLEDA database. The rest of the plots include all the sample. 
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(Left) Contours of integrated HI intensity (zeroth moment) overlaid on the DSS2 image, and (Right) map of intensity-weighted HI mean 
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left corner of the lower panels is the half-power width of the synthesized beam, and has a size of 58^' x 53'^ (full resolution) and 126'' x 113'' 
(the same spatial resolution as the Seyfert sample). 



14 



Tang et al. 



o 
o 
o 

CM 
-3 



o 

< 



o 

UJ 
Q 



36 20 
18 
16 
14 
12 
10 

36 20 
18 
16 
14 
12 
10 

36 20 
18 
16 
14 
12 
10 

36 20 
18 
16 
14 
12 
10 

36 20 
18 
16 
14 
12 
10 



2673.3 KM/S. 

: + : 

■0 


2652.3 KM/S. 

: + : 

- C\ 


2631.4 KM/S. 

o 


2610.4 KM/S. 


2589.4 KM/S. 


2568.4 KM/S. 


■1 h- 

2547.5J<M/S. 

1 1 1 


2526.5 KM/S. 


2505.5 KM/S. 


2484.6 KM/S. 

■ w ■ 


1 

2 


!463.6 KM/S. 


2442.7 KM/S. 

/^)) 


2421 .7 KM/S. 


2400.8 KM/S. 


2379.8 KM/S. 


.' 2358.9 KWi/S. 


2338.0 KM/S. 


2317.0 KM/S. 


2296.1 KM/S. 


2275.2 KM/S. 

: + : 



0813 30 



00 12 30 



0813 30 00 12 30 

RIGHT ASCENSION (J2000) 



Fig. 3. — HI channel maps of NGC 2543 with spatial resolution the same as the Seyfert sample. Contour levels are plotted at —3, 3, 6, 
9, 12, 15, 18 X 0.18 mjy beam~^ (l*^)? which corresponds to a HI column density of 2.0 xlO^^ cm~^. The central heliocentric velocity is 
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Fig. 6. — HI channel maps of NGC 2551 with spatial resolution the same as the Seyfert sample. Contour levels are plotted at —3, 3, 
6, 9, 12, 15, 18 X 0.18 mjy beam~^ (l*^)? which corresponds to a HI column density of 2.0 xlO^^ cm~^. The cross marks the position of 
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Fig. 7. — Upper panel: Optical images of NGC 3094 (Control sample) from DSS2. Middle panels: (Left) Contours of zeroth moment 
overlaid on the DSS2 image and (Right) first moment map with full resolution. Lower panels: (Left) Contours of zeroth moment and (Right) 
first moment map with the same spatial resolution as the Seyfert sample. In the zeroth moment maps with full resolution, contours are 
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the Seyfert sample, contours are plotted at 3, 20, 40, 60, 80 x 71.7 mJy beam~^ km s~^(8.5 x 10^^ cm~^). The half-power width of the 
synthesized beam has a size of 59'' x 52" (full resolution) and 123'' x 111" (the same spatial resolution as the Seyfert sample). 
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Fig. 8. — HI channel maps of NGC 3094 with spatial resolution the same as the Seyfert sample. Contour levels are plotted at —3, —2, 
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Fig. 9. — Upper panel: Optical images of NGC 3976 (Control sample) from DSS2. Middle panels: (Left) Contours of zeroth moment 
overlaid on the DSS2 image and (Right) first moment map with full resolution. Lower panels: (Left) Contours of zeroth moment and 
(Right) first moment map with the same spatial resolution as the Seyfert sample. In the zeroth moment maps with full resolution, contours 
are plotted at 3, 20, 40, 60, 80 x 28.2 mJy beam"-*^ km s~-'^(l.l x lO-*^^ cm~^). In the zeroth moment maps with the same spatial resolution 
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Fig. 10. — Upper panel: Optical image of IC 2461 (Control sample) from the DSS2. Lower panels: (Left) Contours of zeroth moment 
overlaid on the DSS2 image, and (Right) first moment map. In the zeroth moment map, contours are plotted at 3, 20, 40, 60, 80 X 
26.8 mjy beam~^ km s~^(l.l X 10^^ cm~-^). The half-power width of the synthesized beam has a size of 72" x 55" . 
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Fig. 11. — Upper panel: Optical image of NGC 3835 (Control sample) from the DSS2. Lower panels: (Left) Contours of zeroth 
moment overlaid on the DSS2 image, and (Right) first moment map. In the zeroth moment map, contours are plotted at 3, 10, 20, 30, 
40 X 31.2 mjy beam~^ km s~^ (1.6 x 10^^ cm~^). The half-power width of the synthesized beam has a size of 63^' x 51". 
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Fig. 12. — Upper panel: Optical image of NGC 4067 (Control sample) from the DSS2. Lower panels: (Left) Contours of zeroth 
moment overlaid on the DSS2 image, and (Right) first moment map. In the zeroth moment map, contours are plotted at 3, 10, 20, 30, 
40 X 31.2 mjy beam~^ km s~^ (1.5 x 10^^ cm~^). The half-power width of the synthesized beam has a size of 62'' x 54". 
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Fig. 13. — Upper panel: Optical image of NGC 4256 (Control sample) from the DSS2. Lower panels: (Left) Contours of zeroth 
moment overlaid on the DSS2 image, and (Right) first moment map. In the zeroth moment map, contours are plotted at 3, 6, 9, 12, 
15 X 38.0 mjy beam~^ km s~^ (1.8 x 10^^ cm~^). The half-power width of the synthesized beam has a size of 67" x 51". 
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Fig. 14. — Upper panel: Optical image of NGC 4275 (Control sample) from the DSS2. Lower panels: (Left) Contours of zeroth 
moment overlaid on the DSS2 image, and (Right) first moment map. In the zeroth moment map, contours are plotted at 3, 10, 20, 30, 
40 X 23.2 mjy beam~^ km s~^ (1.1 X 10^^ cm~^). The half-power width of the synthesized beam has a size of 63^' x 55". 
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Fig. 15. — Upper panel: Optical image of NGC 4351 (Control sample) from the DSS2. Lower panels: (Left) Contours of zeroth 
moment overlaid on the DSS2 image, and (Right) first moment map. In the zeroth moment map, contours are plotted at 3, 10, 20, 30, 
40 X 20.8 mjy beam~^ km s~^ (1.1 X 10^^ cm~^). The half-power width of the synthesized beam has a size of 59'' x 54". 
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Fig. 16. — Upper panel: Optical image of NGC 4384 (Control sample) from the DSS2. Lower panels: (Left) Contours of zeroth 
moment overlaid on the DSS2 image, and (Right) first moment map. In the zeroth moment map, contours are plotted at 3, 10, 20, 30, 
40 X 18.0 mjy beam~^ km s~^ (7.8 x 10^^ cm"^). The half-power width of the synthesized beam has a size of 66^^ x 56^^ 
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Fig. 17. — Upper panel: Optical image of NGC 4470 (Control sample) from the DSS2. Lower panels: (Left) Contours of zeroth 
moment overlaid on the DSS2 image, and (Right) first moment map. In the zeroth moment map, contours are plotted at 3, 20, 40, 60, 
80 X 17.0 mjy beam~^ km s~^ (8.4 x 10^^ cm"^). The half-power width of the synthesized beam has a size of 62" x 53^^ 
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Fig. 18. — Upper panel: Optical image of NGC 4513 (Control sample) from the DSS2. Lower panels: (Left) Contours of zeroth 
moment overlaid on the DSS2 image, and (Right) first moment map. In the zeroth moment map, contours are plotted at 3, 6, 9, 12, 
15 X 31.2 mjy beam~^ km s~^ (1.2 x 10^^ cm"^). The half-power width of the synthesized beam has a size of 73" x 57". 
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Fig. 19. — Upper panel: Optical image of NGC 4567 (Control sample) from the DSS2. Lower panels: (Left) Contours of zeroth 
moment overlaid on the DSS2 image, and (Right) first moment map. In the zeroth moment map, contours are plotted at 3, 20, 40, 60, 
80 X 24.8 mjy beam~^ km s~^ (1.1 x 10^^ cm~^). The half-power width of the synthesized beam has a size of 64^^ x 57". 
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Fig. 20. — Upper panel: Optical image of NGC 4591 (Control sample) from the DSS2. Lower panels: (Left) Contours of zeroth 
moment overlaid on the DSS2 image, and (Right) first moment map. In the zeroth moment map, contours are plotted at 3, 10, 20, 30, 
40 X 24.0 mjy beam~^ km s~^ (1.1 X 10^^ cm~^). The half-power width of the synthesized beam has a size of 67" x 54". 
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Fig. 21. — Upper panel: Optical image of NGC 4814 (Control sample) from the DSS2. Lower panels: (Left) Contours of zeroth 
moment overlaid on the DSS2 image, and (Right) first moment map. In the zeroth moment map, contours are plotted at 3, 10, 20, 30, 
40 X 46.1 mjy beam~^ km s~^ (2.0 x 10^^ cm~^). The half-power width of the synthesized beam has a size of 71" x 54". 
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Fig. 22. — Upper panel: Optical image of NGC 4964 (Control sample) from the DSS2. Lower panels: (Left) Contours of zeroth 
moment overlaid on the DSS2 image, and (Right) first moment map. In the zeroth moment map, contours are plotted at 3, 10, 20, 30, 
40 X 29.4 mjy beam~^ km s~^ (1.4 x 10^^ cm"^). The half-power width of the synthesized beam has a size of 67" x 52^^ 
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Fig. 23. — Upper panel: Optical image of NGC 5289 (Control sample) and NGC 5290 from the DSS2. Middle panels: Contours of zeroth 
moment overlaid on the DSS2 image. Lower panels: (Left) Contours of zeroth moment overlaid on the DSS2 image and (Right) first moment 
map with larger field. In the zeroth moment maps, contours are plotted at 3, 10, 20, 30, 40 X 20.8mJybeam~^kms~^ (9.6 X lO^^cm""^). 
The half-power width of the synthesized beam has a size of 62" x 56" . 
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Fig. 24. — Upper panel: Optical image of NGC 5350 (Control sample) from the DSS2. Lower panels: (Left) Contours of zeroth 
moment overlaid on the DSS2 image, and (Right) first moment map. In the zeroth moment map, contours are plotted at 3, 10, 20, 30, 
40 X 31.8 mjy beam~^ km s~^ (1.5 x 10^^ cm~-^). The half-power width of the synthesized beam has a size of 61" x 56" . 
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Fig. 25. — Upper panel: Optical image of NGC 5375 (Control sample) from the DSS2. Lower panels: (Left) Contours of zeroth 
moment overlaid on the DSS2 image, and (Right) first moment map. In the zeroth moment map, contours are plotted at 3, 10, 20, 30, 
40 X 32.3 mjy beam~^ km s~^ (1.8 x 10^^ cm"^). The half-power width of the synthesized beam has a size of 57" x 52^^ 
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Fig. 26. — Distribution in HI gas masses of ensemble active galaxy sample o f|Kuo et al.| ( [2008 t (solid histogram) and our ensemble control 
sample (hatched histogram). The combined HI gas masses of NGC 4567/NG(J 4568 in the control sample is plotted. 
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Fig. 27. — Cumulative fraction of Seyfert galaxies (solid line) and control sample (dashed line) with interacting neighboring galaxies 
plotted as a function of their projected separations. 



